Ex vivo tissue perturbations coupled to single-cell RNA-seq reveal multilineage cell circuit dynamics in human lung fibrogenesis.

Pulmonary fibrosis develops as a consequence of failed regeneration after injury. Analyzing mechanisms of regeneration and fibrogenesis directly in human tissue has been hampered by the lack of organotypic models and analytical techniques. In this work, we coupled ex vivo cytokine and drug perturbations of human precision-cut lung slices (hPCLS) with single-cell RNA sequencing and induced a multilineage circuit of fibrogenic cell states in hPCLS. We showed that these cell states were highly similar to the in vivo cell circuit in a multicohort lung cell atlas from patients with pulmonary fibrosis. Using micro-CT-staged patient tissues, we characterized the appearance and interaction of myofibroblasts, an ectopic endothelial cell state, and basaloid epithelial cells in the thickened alveolar septum of early-stage lung fibrosis. Induction of these states in the hPCLS model provided evidence that the basaloid cell state was derived from alveolar type 2 cells, whereas the ectopic endothelial cell state emerged from capillary cell plasticity. Cell-cell communication routes in patients were largely conserved in hPCLS, and antifibrotic drug treatments showed highly cell type-specific effects. Our work provides an experimental framework for perturbational single-cell genomics directly in human lung tissue that enables analysis of tissue homeostasis, regeneration, and pathology. We further demonstrate that hPCLS offer an avenue for scalable, high-resolution drug testing to accelerate antifibrotic drug development and translation.

An integrated cell atlas of the lung in health and disease.

Single-cell technologies have transformed our understanding of human tissues. Yet, studies typically capture only a limited number of donors and disagree on cell type definitions. Integrating many single-cell datasets can address these limitations of individual studies and capture the variability present in the population. Here we present the integrated Human Lung Cell Atlas (HLCA), combining 49 datasets of the human respiratory system into a single atlas spanning over 2.4 million cells from 486 individuals. The HLCA presents a consensus cell type re-annotation with matching marker genes, including annotations of rare and previously undescribed cell types. Leveraging the number and diversity of individuals in the HLCA, we identify gene modules that are associated with demographic covariates such as age, sex and body mass index, as well as gene modules changing expression along the proximal-to-distal axis of the bronchial tree. Mapping new data to the HLCA enables rapid data annotation and interpretation. Using the HLCA as a reference for the study of disease, we identify shared cell states across multiple lung diseases, including SPP1+ profibrotic monocyte-derived macrophages in COVID-19, pulmonary fibrosis and lung carcinoma. Overall, the HLCA serves as an example for the development and use of large-scale, cross-dataset organ atlases within the Human Cell Atlas.

The role of mtDNA haplogroups on metabolic features in narcolepsy type 1.

Narcolepsy type 1 (NT1) is due to selective loss of hypocretin (hcrt)-producing-neurons. Hcrt is a neuropeptide regulating the sleep/wake cycle, as well as feeding behavior. A subset of NT1 patients become overweight/obese, with a dysmetabolic phenotype. We hypothesized that mitochondrial DNA (mtDNA) sequence variation might contribute to the metabolic features in NT1 and we undertook an exploratory survey of mtDNA haplogroups in a cohort of well-characterized patients. We studied 246 NT1 Italian patients, fully defined for their metabolic features, including obesity, hypertension, low HDL, hypertriglyceridemia and hyperglycemia. For haplogroup assignment, the mtDNA control region was sequenced in combination with an assessment of diagnostic markers in the coding region. NT1 patients displayed the same mtDNA haplogroups (H, HV, J, K, T, U) frequency as those reported in the general Italian population. The majority of NT1 patients (64%) were overweight: amongst these, 35% were obese, 48% had low HDL cholesterol levels, and 31% had hypertriglyceridemia. We identified an association between haplogroups J, K and hypertriglyceridemia (P = 0.03, 61.5% and 61.5%, respectively vs. 31.3% of the whole sample) and after correction for age and sex, we observed a reduction of these associations (OR = 3.65, 95%CI = 0.76-17.5, p = 0.106 and 1.73, 0.52-5.69, p = 0.368, respectively). The low HDL level showed a trend for association with haplogroup J (P = 0.09, 83.3% vs. 47.4% of the whole sample) and after correction we observed an OR = 6.73, 95%CI = 0.65-69.9, p = 0.110. Our study provides the first indication that mtDNA haplogroups J and K can modulate metabolic features of NT1 patients, linking mtDNA variation to the dysmetabolic phenotype in NT1.

The transcription factors Runx3 and ThPOK cross-regulate acquisition of cytotoxic function by human Th1 lymphocytes.

Cytotoxic CD4 (CD4CTX) T cells are emerging as an important component of antiviral and antitumor immunity, but the molecular basis of their development remains poorly understood. In the context of human cytomegalovirus infection, a significant proportion of CD4 T cells displays cytotoxic functions. We observed that the transcriptional program of these cells was enriched in CD8 T cell lineage genes despite the absence of ThPOK downregulation. We further show that establishment of CD4CTX-specific transcriptional and epigenetic programs occurred in a stepwise fashion along the Th1-differentiation pathway. In vitro, prolonged activation of naive CD4 T cells in presence of Th1 polarizing cytokines led to the acquisition of perforin-dependent cytotoxic activity. This process was dependent on the Th1 transcription factor Runx3 and was limited by the sustained expression of ThPOK. This work elucidates the molecular program of human CD4CTX T cells and identifies potential targets for immunotherapy against viral infections and cancer.

Origin and spread of human mitochondrial DNA haplogroup U7.

Human mitochondrial DNA haplogroup U is among the initial maternal founders in Southwest Asia and Europe and one that best indicates matrilineal genetic continuity between late Pleistocene hunter-gatherer groups and present-day populations of Europe. While most haplogroup U subclades are older than 30 thousand years, the comparatively recent coalescence time of the extant variation of haplogroup U7 (~16-19 thousand years ago) suggests that its current distribution is the consequence of more recent dispersal events, despite its wide geographical range across Europe, the Near East and South Asia. Here we report 267 new U7 mitogenomes that - analysed alongside 100 published ones - enable us to discern at least two distinct temporal phases of dispersal, both of which most likely emanated from the Near East. The earlier one began prior to the Holocene (~11.5 thousand years ago) towards South Asia, while the later dispersal took place more recently towards Mediterranean Europe during the Neolithic (~8 thousand years ago). These findings imply that the carriers of haplogroup U7 spread to South Asia and Europe before the suggested Bronze Age expansion of Indo-European languages from the Pontic-Caspian Steppe region.

Mitogenomes from two uncommon haplogroups mark late glacial/postglacial expansions from the near east and neolithic dispersals within Europe.

The current human mitochondrial (mtDNA) phylogeny does not equally represent all human populations but is biased in favour of representatives originally from north and central Europe. This especially affects the phylogeny of some uncommon West Eurasian haplogroups, including I and W, whose southern European and Near Eastern components are very poorly represented, suggesting that extensive hidden phylogenetic substructure remains to be uncovered. This study expanded and re-analysed the available datasets of I and W complete mtDNA genomes, reaching a comprehensive 419 mitogenomes, and searched for precise correlations between the ages and geographical distributions of their numerous newly identified subclades with events of human dispersal which contributed to the genetic formation of modern Europeans. Our results showed that haplogroups I (within N1a1b) and W originated in the Near East during the Last Glacial Maximum or pre-warming period (the period of gradual warming between the end of the LGM, ∼19 ky ago, and the beginning of the first main warming phase, ∼15 ky ago) and, like the much more common haplogroups J and T, may have been involved in Late Glacial expansions starting from the Near East. Thus our data contribute to a better definition of the Late and postglacial re-peopling of Europe, providing further evidence for the scenario that major population expansions started after the Last Glacial Maximum but before Neolithic times, but also evidencing traces of diffusion events in several I and W subclades dating to the European Neolithic and restricted to Europe.

Reconciling migration models to the Americas with the variation of North American native mitogenomes.

In this study we evaluated migration models to the Americas by using the information contained in native mitochondrial genomes (mitogenomes) from North America. Molecular and phylogeographic analyses of B2a mitogenomes, which are absent in Eskimo-Aleut and northern Na-Dene speakers, revealed that this haplogroup arose in North America ∼11-13 ka from one of the founder Paleo-Indian B2 mitogenomes. In contrast, haplogroup A2a, which is typical of Eskimo-Aleuts and Na-Dene, but also present in the easternmost Siberian groups, originated only 4-7 ka in Alaska, led to the first Paleo-Eskimo settlement of northern Canada and Greenland, and contributed to the formation of the Na-Dene gene pool. However, mitogenomes also show that Amerindians from northern North America, without any distinction between Na-Dene and non-Na-Dene, were heavily affected by an additional and distinctive Beringian genetic input. In conclusion, most mtDNA variation (along the double-continent) stems from the first wave from Beringia, which followed the Pacific coastal route. This was accompanied or followed by a second inland migratory event, marked by haplogroups X2a and C4c, which affected all Amerindian groups of Northern North America. Much later, the ancestral A2a carriers spread from Alaska, undertaking both a westward migration to Asia and an eastward expansion into the circumpolar regions of Canada. Thus, the first American founders left the greatest genetic mark but the original maternal makeup of North American Natives was subsequently reshaped by additional streams of gene flow and local population dynamics, making a three-wave view too simplistic.

Uniparental genetic heritage of belarusians: encounter of rare middle eastern matrilineages with a central European mitochondrial DNA pool.

Ethnic Belarusians make up more than 80% of the nine and half million people inhabiting the Republic of Belarus. Belarusians together with Ukrainians and Russians represent the East Slavic linguistic group, largest both in numbers and territory, inhabiting East Europe alongside Baltic-, Finno-Permic- and Turkic-speaking people. Till date, only a limited number of low resolution genetic studies have been performed on this population. Therefore, with the phylogeographic analysis of 565 Y-chromosomes and 267 mitochondrial DNAs from six well covered geographic sub-regions of Belarus we strove to complement the existing genetic profile of eastern Europeans. Our results reveal that around 80% of the paternal Belarusian gene pool is composed of R1a, I2a and N1c Y-chromosome haplogroups - a profile which is very similar to the two other eastern European populations - Ukrainians and Russians. The maternal Belarusian gene pool encompasses a full range of West Eurasian haplogroups and agrees well with the genetic structure of central-east European populations. Our data attest that latitudinal gradients characterize the variation of the uniparentally transmitted gene pools of modern Belarusians. In particular, the Y-chromosome reflects movements of people in central-east Europe, starting probably as early as the beginning of the Holocene. Furthermore, the matrilineal legacy of Belarusians retains two rare mitochondrial DNA haplogroups, N1a3 and N3, whose phylogeographies were explored in detail after de novo sequencing of 20 and 13 complete mitogenomes, respectively, from all over Eurasia. Our phylogeographic analyses reveal that two mitochondrial DNA lineages, N3 and N1a3, both of Middle Eastern origin, might mark distinct events of matrilineal gene flow to Europe: during the mid-Holocene period and around the Pleistocene-Holocene transition, respectively.

Autosomal and uniparental portraits of the native populations of Sakha (Yakutia): implications for the peopling of Northeast Eurasia.

BACKGROUND: Sakha--an area connecting South and Northeast Siberia--is significant for understanding the history of peopling of Northeast Eurasia and the Americas. Previous studies have shown a genetic contiguity between Siberia and East Asia and the key role of South Siberia in the colonization of Siberia. RESULTS: We report the results of a high-resolution phylogenetic analysis of 701 mtDNAs and 318 Y chromosomes from five native populations of Sakha (Yakuts, Evenks, Evens, Yukaghirs and Dolgans) and of the analysis of more than 500,000 autosomal SNPs of 758 individuals from 55 populations, including 40 previously unpublished samples from Siberia. Phylogenetically terminal clades of East Asian mtDNA haplogroups C and D and Y-chromosome haplogroups N1c, N1b and C3, constituting the core of the gene pool of the native populations from Sakha, connect Sakha and South Siberia. Analysis of autosomal SNP data confirms the genetic continuity between Sakha and South Siberia. Maternal lineages D5a2a2, C4a1c, C4a2, C5b1b and the Yakut-specific STR sub-clade of Y-chromosome haplogroup N1c can be linked to a migration of Yakut ancestors, while the paternal lineage C3c was most likely carried to Sakha by the expansion of the Tungusic people. MtDNA haplogroups Z1a1b and Z1a3, present in Yukaghirs, Evens and Dolgans, show traces of different and probably more ancient migration(s). Analysis of both haploid loci and autosomal SNP data revealed only minor genetic components shared between Sakha and the extreme Northeast Siberia. Although the major part of West Eurasian maternal and paternal lineages in Sakha could originate from recent admixture with East Europeans, mtDNA haplogroups H8, H20a and HV1a1a, as well as Y-chromosome haplogroup J, more probably reflect an ancient gene flow from West Eurasia through Central Asia and South Siberia. CONCLUSIONS: Our high-resolution phylogenetic dissection of mtDNA and Y-chromosome haplogroups as well as analysis of autosomal SNP data suggests that Sakha was colonized by repeated expansions from South Siberia with minor gene flow from the Lower Amur/Southern Okhotsk region and/or Kamchatka. The minor West Eurasian component in Sakha attests to both recent and ongoing admixture with East Europeans and an ancient gene flow from West Eurasia.

Small effective population size and genetic homogeneity in the Val Borbera isolate.

Population isolates are a valuable resource for medical genetics because of their reduced genetic, phenotypic and environmental heterogeneity. Further, extended linkage disequilibrium (LD) allows accurate haplotyping and imputation. In this study, we use nuclear and mitochondrial DNA data to determine to what extent the geographically isolated population of the Val Borbera valley also presents features of genetic isolation. We performed a comparative analysis of population structure and estimated effective population size exploiting LD data. We also evaluated haplotype sharing through the analysis of segments of autozygosity. Our findings reveal that the valley has features characteristic of a genetic isolate, including reduced genetic heterogeneity and reduced effective population size. We show that this population has been subject to prolonged genetic drift and thus we expect many variants that are rare in the general population to reach significant frequency values in the valley, making this population suitable for the identification of rare variants underlying complex traits.

Rare primary mitochondrial DNA mutations and probable synergistic variants in Leber's hereditary optic neuropathy.

BACKGROUND: Leber's hereditary optic neuropathy (LHON) is a maternally inherited blinding disorder, which in over 90% of cases is due to one of three primary mitochondrial DNA (mtDNA) point mutations (m.11778G>A, m.3460G>A and m.14484T>C, respectively in MT-ND4, MT-ND1 and MT-ND6 genes). However, the spectrum of mtDNA mutations causing the remaining 10% of cases is only partially and often poorly defined. METHODOLOGY/PRINCIPAL FINDINGS: In order to improve such a list of pathological variants, we completely sequenced the mitochondrial genomes of suspected LHON patients from Italy, France and Germany, lacking the three primary common mutations. Phylogenetic and conservation analyses were performed. Sixteen mitochondrial genomes were found to harbor at least one of the following nine rare LHON pathogenic mutations in genes MT-ND1 (m.3700G>A/p.A132T, m.3733G>A-C/p.E143K-Q, m.4171C>A/p.L289M), MT-ND4L (m.10663T>C/p.V65A) and MT-ND6 (m.14459G>A/p.A72V, m.14495A>G/p.M64I, m.14482C>A/p.L60S, and m.14568C>T/p.G36S). Phylogenetic analyses revealed that these substitutions were due to independent events on different haplogroups, whereas interspecies comparisons showed that they affected conserved amino acid residues or domains in the ND subunit genes of complex I. CONCLUSIONS/SIGNIFICANCE: Our findings indicate that these nine substitutions are all primary LHON mutations. Therefore, despite their relative low frequency, they should be routinely tested for in all LHON patients lacking the three common mutations. Moreover, our sequence analysis confirms the major role of haplogroups J1c and J2b (over 35% in our probands versus 6% in the general population of Western Europe) and other putative synergistic mtDNA variants in LHON expression.

Ancient migratory events in the Middle East: new clues from the Y-chromosome variation of modern Iranians.

Knowledge of high resolution Y-chromosome haplogroup diversification within Iran provides important geographic context regarding the spread and compartmentalization of male lineages in the Middle East and southwestern Asia. At present, the Iranian population is characterized by an extraordinary mix of different ethnic groups speaking a variety of Indo-Iranian, Semitic and Turkic languages. Despite these features, only few studies have investigated the multiethnic components of the Iranian gene pool. In this survey 938 Iranian male DNAs belonging to 15 ethnic groups from 14 Iranian provinces were analyzed for 84 Y-chromosome biallelic markers and 10 STRs. The results show an autochthonous but non-homogeneous ancient background mainly composed by J2a sub-clades with different external contributions. The phylogeography of the main haplogroups allowed identifying post-glacial and Neolithic expansions toward western Eurasia but also recent movements towards the Iranian region from western Eurasia (R1b-L23), Central Asia (Q-M25), Asia Minor (J2a-M92) and southern Mesopotamia (J1-Page08). In spite of the presence of important geographic barriers (Zagros and Alborz mountain ranges, and the Dasht-e Kavir and Dash-e Lut deserts) which may have limited gene flow, AMOVA analysis revealed that language, in addition to geography, has played an important role in shaping the nowadays Iranian gene pool. Overall, this study provides a portrait of the Y-chromosomal variation in Iran, useful for depicting a more comprehensive history of the peoples of this area as well as for reconstructing ancient migration routes. In addition, our results evidence the important role of the Iranian plateau as source and recipient of gene flow between culturally and genetically distinct populations.

Distinguishing the co-ancestries of haplogroup G Y-chromosomes in the populations of Europe and the Caucasus.

Haplogroup G, together with J2 clades, has been associated with the spread of agriculture, especially in the European context. However, interpretations based on simple haplogroup frequency clines do not recognize underlying patterns of genetic diversification. Although progress has been recently made in resolving the haplogroup G phylogeny, a comprehensive survey of the geographic distribution patterns of the significant sub-clades of this haplogroup has not been conducted yet. Here we present the haplogroup frequency distribution and STR variation of 16 informative G sub-clades by evaluating 1472 haplogroup G chromosomes belonging to 98 populations ranging from Europe to Pakistan. Although no basal G-M201* chromosomes were detected in our data set, the homeland of this haplogroup has been estimated to be somewhere nearby eastern Anatolia, Armenia or western Iran, the only areas characterized by the co-presence of deep basal branches as well as the occurrence of high sub-haplogroup diversity. The P303 SNP defines the most frequent and widespread G sub-haplogroup. However, its sub-clades have more localized distribution with the U1-defined branch largely restricted to Near/Middle Eastern and the Caucasus, whereas L497 lineages essentially occur in Europe where they likely originated. In contrast, the only U1 representative in Europe is the G-M527 lineage whose distribution pattern is consistent with regions of Greek colonization. No clinal patterns were detected suggesting that the distributions are rather indicative of isolation by distance and demographic complexities.

Mitochondrial DNA signals of late glacial recolonization of Europe from near eastern refugia.

Human populations, along with those of many other species, are thought to have contracted into a number of refuge areas at the height of the last Ice Age. European populations are believed to be, to a large extent, the descendants of the inhabitants of these refugia, and some extant mtDNA lineages can be traced to refugia in Franco-Cantabria (haplogroups H1, H3, V, and U5b1), the Italian Peninsula (U5b3), and the East European Plain (U4 and U5a). Parts of the Near East, such as the Levant, were also continuously inhabited throughout the Last Glacial Maximum, but unlike western and eastern Europe, no archaeological or genetic evidence for Late Glacial expansions into Europe from the Near East has hitherto been discovered. Here we report, on the basis of an enlarged whole-genome mitochondrial database, that a substantial, perhaps predominant, signal from mitochondrial haplogroups J and T, previously thought to have spread primarily from the Near East into Europe with the Neolithic population, may in fact reflect dispersals during the Late Glacial period, ∼19-12 thousand years (ka) ago.

Mitochondrial genomes from modern horses reveal the major haplogroups that underwent domestication.

Archaeological and genetic evidence concerning the time and mode of wild horse (Equus ferus) domestication is still debated. High levels of genetic diversity in horse mtDNA have been detected when analyzing the control region; recurrent mutations, however, tend to blur the structure of the phylogenetic tree. Here, we brought the horse mtDNA phylogeny to the highest level of molecular resolution by analyzing 83 mitochondrial genomes from modern horses across Asia, Europe, the Middle East, and the Americas. Our data reveal 18 major haplogroups (A-R) with radiation times that are mostly confined to the Neolithic and later periods and place the root of the phylogeny corresponding to the Ancestral Mare Mitogenome at ~130-160 thousand years ago. All haplogroups were detected in modern horses from Asia, but F was only found in E. przewalskii--the only remaining wild horse. Therefore, a wide range of matrilineal lineages from the extinct E. ferus underwent domestication in the Eurasian steppes during the Eneolithic period and were transmitted to modern E. caballus breeds. Importantly, now that the major horse haplogroups have been defined, each with diagnostic mutational motifs (in both the coding and control regions), these haplotypes could be easily used to (i) classify well-preserved ancient remains, (ii) (re)assess the haplogroup variation of modern breeds, including Thoroughbreds, and (iii) evaluate the possible role of mtDNA backgrounds in racehorse performance.

Bulgarians vs the other European populations: a mitochondrial DNA perspective.

To define the matrilineal relationships between Bulgarians and other European populations, we have evaluated the mitochondrial DNA (mtDNA) variation in a sample of 855 Bulgarian subjects from the mtDNA perspective. The molecular survey was performed by sequencing ∼750 bp of the control region, which resulted in 557 different haplotypes, and by a subsequent restriction fragment length polymorphism analysis to confirm haplogroup/subhaplogroup affiliation. The classification was carried out according to the most updated criteria as reported by van Oven and Kayser (Hum Mutat 30:386-394, 2009), allowing the identification of 45 mitochondrial clades. The observed pattern of mtDNA variation indicates that the Bulgarian mitochondrial pool is geographically homogeneous across the country, and that is characterized by an overall extremely high frequency of western Eurasian lineages. In the principal component analysis, Bulgarians locate in an intermediate position between Eastern European and Mediterranean populations, which is in agreement with historical events. Thus, while the Mediterranean legacy could be attributed to the Thracians, indigenous people that firstly inhabited the Balkans, the Eastern contribution is likely due to the Proto-Bulgarians originating from the Middle East and to the Slavs migrating from northeast Europe.

Mitochondrial haplogroup C4c: a rare lineage entering America through the ice-free corridor?

Recent analyses of mitochondrial genomes from Native Americans have brought the overall number of recognized maternal founding lineages from just four to a current count of 15. However, because of their relative low frequency, almost nothing is known for some of these lineages. This leaves a considerable void in understanding the events that led to the colonization of the Americas following the Last Glacial Maximum (LGM). In this study, we identified and completely sequenced 14 mitochondrial DNAs belonging to one extremely rare Native American lineage known as haplogroup C4c. Its age and geographical distribution raise the possibility that C4c marked the Paleo-Indian group(s) that entered North America from Beringia through the ice-free corridor between the Laurentide and Cordilleran ice sheets. The similarities in ages andgeographical distributions for C4c and the previously analyzed X2a lineage provide support to the scenario of a dual origin for Paleo-Indians. Taking into account that C4c is deeply rooted in the Asian portion of the mtDNA phylogeny and is indubitably of Asian origin, the finding that C4c and X2a are characterized by parallel genetic histories definitively dismisses the controversial hypothesis of an Atlantic glacial entry route into North America.

In search of the genetic footprints of Sumerians: a survey of Y-chromosome and mtDNA variation in the Marsh Arabs of Iraq.

BACKGROUND: For millennia, the southern part of the Mesopotamia has been a wetland region generated by the Tigris and Euphrates rivers before flowing into the Gulf. This area has been occupied by human communities since ancient times and the present-day inhabitants, the Marsh Arabs, are considered the population with the strongest link to ancient Sumerians. Popular tradition, however, considers the Marsh Arabs as a foreign group, of unknown origin, which arrived in the marshlands when the rearing of water buffalo was introduced to the region. RESULTS: To shed some light on the paternal and maternal origin of this population, Y chromosome and mitochondrial DNA (mtDNA) variation was surveyed in 143 Marsh Arabs and in a large sample of Iraqi controls. Analyses of the haplogroups and sub-haplogroups observed in the Marsh Arabs revealed a prevalent autochthonous Middle Eastern component for both male and female gene pools, with weak South-West Asian and African contributions, more evident in mtDNA. A higher male than female homogeneity is characteristic of the Marsh Arab gene pool, likely due to a strong male genetic drift determined by socio-cultural factors (patrilocality, polygamy, unequal male and female migration rates). CONCLUSIONS: Evidence of genetic stratification ascribable to the Sumerian development was provided by the Y-chromosome data where the J1-Page08 branch reveals a local expansion, almost contemporary with the Sumerian City State period that characterized Southern Mesopotamia. On the other hand, a more ancient background shared with Northern Mesopotamia is revealed by the less represented Y-chromosome lineage J1-M267*. Overall our results indicate that the introduction of water buffalo breeding and rice farming, most likely from the Indian sub-continent, only marginally affected the gene pool of autochthonous people of the region. Furthermore, a prevalent Middle Eastern ancestry of the modern population of the marshes of southern Iraq implies that if the Marsh Arabs are descendants of the ancient Sumerians, also the Sumerians were most likely autochthonous and not of Indian or South Asian ancestry.

Mitochondrial DNA backgrounds might modulate diabetes complications rather than T2DM as a whole.

Mitochondrial dysfunction has been implicated in rare and common forms of type 2 diabetes (T2DM). Additionally, rare mitochondrial DNA (mtDNA) mutations have been shown to be causal for T2DM pathogenesis. So far, many studies have investigated the possibility that mtDNA variation might affect the risk of T2DM, however, when found, haplogroup association has been rarely replicated, even in related populations, possibly due to an inadequate level of haplogroup resolution. Effects of mtDNA variation on diabetes complications have also been proposed. However, additional studies evaluating the mitochondrial role on both T2DM and related complications are badly needed. To test the hypothesis of a mitochondrial genome effect on diabetes and its complications, we genotyped the mtDNAs of 466 T2DM patients and 438 controls from a regional population of central Italy (Marche). Based on the most updated mtDNA phylogeny, all 904 samples were classified into 57 different mitochondrial sub-haplogroups, thus reaching an unprecedented level of resolution. We then evaluated whether the susceptibility of developing T2DM or its complications differed among the identified haplogroups, considering also the potential effects of phenotypical and clinical variables. MtDNA backgrounds, even when based on a refined haplogroup classification, do not appear to play a role in developing T2DM despite a possible protective effect for the common European haplogroup H1, which harbors the G3010A transition in the MTRNR2 gene. In contrast, our data indicate that different mitochondrial haplogroups are significantly associated with an increased risk of specific diabetes complications: H (the most frequent European haplogroup) with retinopathy, H3 with neuropathy, U3 with nephropathy, and V with renal failure.

Mitochondrial haplogroup H1 in north Africa: an early holocene arrival from Iberia.

The Tuareg of the Fezzan region (Libya) are characterized by an extremely high frequency (61%) of haplogroup H1, a mitochondrial DNA (mtDNA) haplogroup that is common in all Western European populations. To define how and when H1 spread from Europe to North Africa up to the Central Sahara, in Fezzan, we investigated the complete mitochondrial genomes of eleven Libyan Tuareg belonging to H1. Coalescence time estimates suggest an arrival of the European H1 mtDNAs at about 8,000-9,000 years ago, while phylogenetic analyses reveal three novel H1 branches, termed H1v, H1w and H1x, which appear to be specific for North African populations, but whose frequencies can be extremely different even in relatively close Tuareg villages. Overall, these findings support the scenario of an arrival of haplogroup H1 in North Africa from Iberia at the beginning of the Holocene, as a consequence of the improvement in climate conditions after the Younger Dryas cold snap, followed by in situ formation of local H1 sub-haplogroups. This process of autochthonous differentiation continues in the Libyan Tuareg who, probably due to isolation and recent founder events, are characterized by village-specific maternal mtDNA lineages.